Smallest Dirac Eigenvalue Distribution from Random Matrix Theory 



OO 

S3 
S3 



(N 
> 

o 
o 
m 
o 

oo 
On 



i 

^ ■ 



X 



Shinsuke M. Nishigaki 1 '*, Poul H. Dam gaard and Tilo Wettig 
1 Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030, USA 

2 The Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
3 ' Institut fur Theoretische Physik, Technische Universitat Miinchen, D-85747 G 'arching, Germany 

(March 2, 1998) 

We derive the hole probability and the distribution of the smallest eigenvalue of chiral hermitian 
random matrices corresponding to Dirac operators coupled to massive quarks in QCD. They are 
expressed in terms of the QCD partition function in the mesoscopic regime. Their universality is 
explicitly related to that of the microscopic massive Bessel kernel. 
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There has long been an attractive idea that the low- 
energy physics of a complex system can be described by 
a simple effective theory which respects the global sym- 
metries of the original system. As an example, the quan- 
tum spectral statistics of a classically chaotic system is 
believed to be described by a random matrix theory be- 
longing to the same universality class as the former |Q. 
One new manifestation of essentially the same idea is the 
recent observation that QCD Dirac operator spectra on 
the scale A = 0(1/P 4 ) (where V" 4 is the space-time vol- 
ume) measured in lattice Monte Carlo simulations || are 
in excellent agreement with the predictions from those 
large- TV random matrix theories that share the same 
global symmetries as QCD. The suitably rescaled (micro- 
scopic) spectral correlation functions thus seem to pro- 
vide exact finite-size scaling functions for QCD in a finite 
volume. Very recently, the microscopic spectral corre- 
lators have been calculated from random matrix theo- 
ries that include the effect of fermion determinants with 
masses m ~ 0(1/14) (see also ||). When A and 

m are measured in units of the mean level spacing at 
zero virtuality, all the random matrix predictions turn 
out to be universal, i.e., insensitive to the details of the 
random matrix potential j|-[3| ■ Although the question of 
whether or not QCD is included in the same universality 
class cannot be answered by demonstrating the existence 
of the wide range of universality within random matrix 
theories, it provides strong support for the former. 

From the field-theoretic point of view || it would be 
most surprising if these observables would not also be 
computable solely within the framework of finite- volume 
generating functionals (partition functions) for the order 
parameter (ipip). If not, large- N random matrix theory, 
which in principle is foreign to the pertinent field theory 
language, would seem to be a new ingredient required 
to describe the observed spectral correlators. It has re- 
cently been shown that a description entirely in terms of 
finite-volume partition functions is indeed also possible 

& 

In order to confirm by numerical simulations that the 
low-lying spectra of QCD Dirac operators can be de- 
scribed alternatively by large- A random matrix theories, 
it is in practice most convenient to measure the distribu- 



tion of the smallest eigenvalue 0] and compare that to 
the random matrix prediction llj7j]. Since the smallest 



eigenvalue distribution largely consists of the first peak 
of the microscopic spectral density (see Fig. 1, £ = A^A), 
we expect it to be universal. 




FIG. 1. Microscopic spectral density (solid line, normal- 
ized as p s (C^ oo) = I/ 71 ") and smallest eigenvalue distribution 
(dotted line) for the quenched chiral unitary ensemble. 

In fact, the proven universality of the massive kernels 
[Bpl of (chiral) unitary ensembles of random matrices 
guarantees the universality of the hole probability E(s), 
i.e., the probability that the interval [— s,s] is free of 
eigenvalues, and of the smallest eigenvalue distribution 
P(s) — —E'(s) because it is related to the kernel via the 
Fredholm determinant formula 1131, 



E(s) = det(l - K). 



(1) 



Here, K is an integral operator whose kernel is the mi- 
croscopic massive Bessel kernel (Eq. (jlj) below or its 
non-chiral counterpart) over the interval [— s, a]. While 
other universal statistical quantities such as the num- 
ber variance E 2 (s) and the spectral rigidity Aa(s) of the 
eigenvalues in the interval [— s, s] are directly related to 
the kernel by integral transforms, the techniques required 
to compute ([j]) from the kernel are rather involved. In 
the case of the universal (massless) sine kernel describ- 
ing the bulk of unitary ensembles, its eigenmodes are 
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given by the prolate spheroidal functions whose eigenval- 
ues are known. However, the massive generalization of 
these eigenmodes is difficult to construct. The alternative 
technique of determining E(s) as a r-function of a sys- 
tem of integrable partial differential equations Jl3] docs 
not simplify the situation either. Consequently, the ex- 
plicit calculation for the massive case has been done so far 
solely for the chiral Gaussian unitary ensemble 0. The 
result contains an explicit factor of exp(— £ 2 /4) which at 
first sight is due to the Gaussian potential. Since these 
results should be universal, the identical factor should 
also arise from a generic random matrix potential. How 
this happens is not immediately obvious and will be ad- 
dressed in the sequel. We shall provide a concise method 
to circumvent the difficulty in explicitly evaluating (0) 
and to efficiently calculate E(s) for the massive chiral 
unitary ensembles with generic potentials. 

We define the partition function of the chiral unitary 
ensemble 



Z^({m}) 



dMe 



-mtv(M) jj det ( M . 



where M is a 2N x 2N block hermitian matrix 



M 



W 
W f 



(2) 



(3) 



dM is the Haar measure of W, and v(M) is an even 
analytic function. The matrix M models the Dirac op- 
erator for SU(iV c > 3) QCD 4 in the Weyl basis, and N 



can be identified with the spacetime volume V4. The 
integer a corresponds to the number of flavors. The 
case of nonzero topological charge v can be treated 
by introducing \v\ massless flavors (3). One then has 
a = Nf + \v\, where Nj is the number of (massive) flavors 
and m Nf+ i, . . . , rriN f+ \ v \ = 0. 

The (unnormalized) probability of having no eigenval- 
ues in the interval [— s, s) is given by 

E (a) (s; {m}) = / dM e~ Ntr v{M) ]J det ( M + im f ) . 

J |ciffcnvalucs|>s * n 



/=i 



(4) 



It is convenient to change the picture from M to an N x N 
positive definite hermitian matrix H = W^W, 



E^(s;{m})= f dH e- NtlV ^Y[det(H + mj) 

J eigenvalues>s 2 j 

= [ dH e- NtrV ^ + ^Y[det(H + s 2 +m 2 f ), (5) 

J ci£ciivalucs>0 £ 



where V(z 2 ) = 2v(z). Hereafter we write 



rn 



(6) 



We can express (H) in terms of an expectation value 
(•••)„,/ with respect to the measure e - NtlV ( H ) x 
n f det(ff + m' 2 ), 



E^ a \s;{m}) 



f °°dH, 



-NtrV(H+s 2 ) 



Uf det (J? + 



r °° dHe~ NtrV ( H ) Uf det(£T + m' f 2 ) 

Ntv(V(H+s 2 )-V(H))\ Z {a) {{m'}). 



m ' 2\ pea 

/ dHe- m * v WT[det(H + m' f 2 ) 
f ) Jo t 



(7) 



In the following we shall show that the two factors in 
Eq. (0) are piecewise universal in the limit ./V — > 00 with 
£ = Ns and p = Nm fixed finite. 

Due to the large- N factorization of macroscopic corre- 
lation functions of U(./V)-invariant operators O, O' Pjj], 



(OC) = (0)(0'} + 



N 2 



(8) 



the first factor of (|t]) is, in the large- N limit, approxi- 
mated by 



e -N 2 (-Lt r (V(H+s 2 )-V{H))) m , 
e -N 2 (s 2 (±trV'(H))+0(s*))_ 



(9) 



Since the fermion determinant does not contribute to 
the macroscopic correlator in the large- N limit, we have 
dropped the suffix m'. 

Now we change the picture back to M, whose macro- 
scopic spectral density is p(z) = ((l/2N)trS(z — M)) = 



R(z 2 )y/a 2 — z 2 , where [—a, a] is the support of the spec- 
trum of M and R{z 2 ) is an analytic function which de- 
pends on the details of v(M) M. We obtain 



1 v'(M)\ f a , . .v'{z) 

In terms ofjo(z), the leading (of order 0(N 2 )) part of the 
action in (H) is written as 



S = j dz p(z)2v(z) — / dzdw p{z)p(w)2P\ii\z — w\. 

J —a J —a 

(11) 

The second term is the exponentiated Vandermonde de- 
terminant. By substituting the large- N saddle-point 
equation 15 SS/6p(z) = 0, i.e., 
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v'(z) - 2 / dw p(w)P = 0. 

J -a Z ~W 

into (|l0|), we have 

f a 1 1 

Q = 2 dzdw p(z)p(w)P- P 

J -a z z - w 



(12) 



(13) 



That is, 



-Ntr(V(H+s 2 )-V(H)) 



e -(^p(o)C) 2 



(15) 



Using the identity l/(z ± ie) = P(l/z) =F iird(z) and tak- 
ing into account that p(z) is even, we finally obtain 



dz dw p(z)p(w)tt5(z)ttS(z — w) = (irp(O)y 



(14) 



which, after the rescaling £ — ► £ / (2tv p(0)) , universally 
reads exp(— C 2 /4). This answers the question raised in 
the introduction how the Gaussian factor arises from a 
generic potential. 

We denote the microscopic limit of the partition func- 
tion Z( Q )({m} = {#}) by Z^{{p}). It is related to the 
microscopic kernel by the master formula |lOj 



A-W(Ci, C2; mi, •• • = C^Cr^n ^(C 2 + m})(CI + M/) 



^" +2 )«l,iC2,Ml,---^a) 
(/*!,..., Ma) 



(16) 



where C is a normalization constant. This formula is reminiscent of the very definition of the partition function and 
the kernel fil| and is valid in the large- N limit, regardless of whether microscopic or macroscopic eigenvalue and mass 
variables are kept finite. 

Using the technique of orthogonal polynomials and rescaling £ — > ( / (2tt p(0 j) , p — > n/(2np(0)), the left-hand side is 
shown to be universally given by Q 



^ (a) (Ci,C 2 ;M) = C#5 



^ det £y(Ci,C2;{/4) 



Ci 2 - CI 11/ v/(Ci 2 + M/ 2 ) (CI + M/ 2 ) , det ^(M) ' 



(17) 



Ay = /<; Bi 



A; 



(< = 1,. 



(-Ci^y-'Jj-iCCi^) (i = a + l,a+2) ' 



where J and / denote the Bessel function of real and 
imaginary argument, respectively. Therefore, after con- 
tinuing £ — + — ifi, we can derive Z^ 2n \{p}) iteratively 
using (|l6|), starting from Z^ = 1. Z^ 2 ™ -1 ) ({/z}) can be 
obtained by decoupling one of the masses in Z^ 2n \{p}) 
by sending it to infinity. In this way, one obtains univer- 
sally 



z [a) m) 



det Aij{{n}) 

l<z,j<a 



det \x 

1<2,j <a 



2(J-1) 



(18) 



which is identical to the finite-volume QCD partition 
function as calculated from the chiral Lagrangian in the 
"mesoscopic" scaling limit ]9|Jl7| . 

Both factors in (0) being universal, the normalized hole 
probability 



^((•M)=e-^ 4 Z(a)({V ^^ 

K ^ m) ^ (Q) (M) 

as well as the smallest eigenvalue distribution J7J 

p (q) (C;M) = -^ (q) (C;M) 

_ S -^ 2 /4 1<»J<Q 

~ 2 6 det Ay (M) 



(19) 



(20) 



with Cij({p}) = pi 1 Ij+i(pi) are universal. As ex- 
plained above, the general case of nonzero topological 
charge v is obtained by introducing \v\ massless flavors 
(a = JV) + M). 

Except for the Gaussian prefactor the expressions (|l9| ) 
and (^0|) are explicitly given in terms of a finite-volume 
field theory partition function in the mesoscopic scaling 
regime. This indicates that also these quantities can be 
derived directly from field theory in the mesoscopic scal- 
ing regime, without the bypass through random matrix 
theory |fi}. 

The hole probabilities for the orthogonal ((3=1) and 
symplectic (f3 = 4) ensembles are also related to the uni- 
tary kernel K by formulae analogous to Eq. (g) (det± 
represents the determinant projected to even or odd or- 
ders of eigenvalues) (L2) : 



E 1 (s) = det + (l - K), 

E 4 (s) = - (det+(l -K) + (det_(l - K) 



(21a) 
(21b) 



which hold regardless of the fermion determinant. Hence, 
the smallest eigenvalue distributions for these ensembles 
are also guaranteed to be universal. This is consistent 
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with the recently proven universality 01| of the micro- 
scopic kernels for the massless chiral orthogonal and sym- 
plectic ensembles |19 via universal relationships between 
them and the massless chiral unitary kernel. Explicit ex- 
pressions for Ei(s) and E^s) are however not known ex- 
cept for the massless chiral Gaussian/Laguerre case [ g0[ . 
To obtain these quantities for the corresponding massive 
ensembles, our method presented here needs modifica- 
tion because the mapping M i— > H involves a nontrivial 
Jacobian, (det H) 13 / 2-1 . This point will be discussed else- 
where. 
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